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This paper considers the problem of locating a single server on a network, relaxing the assumption 
that the server is always available for service, and explicitly accounting for queueing. The 
resulting queueing-location model allows for an arbitrary number of priority classes. Properties 
of the objective function are developed and algorithms presented for obtaining the optimal 
location on tree and cyclic networks. Sensitivity analysis with respect to the average arrival rate 
of calls is investigated. A numerical example is presented to illustrate the results of this paper. 
The major conclusions of the paper include: (a) the optimal location need not be at a node of 
the network, (b) the optimal location changes as a function of the arrival rate of calls into the 
system, (c) the optimal location is usually different from that obtained by grouping all calls into 
one priority class. 


INTRODUCTION 


In recent papers Berman, Larson, and Chiu [3] and Chiu, Berman, and Larson [5] 
have proposed and analyzed a queueing-location model, which they refer to as the 
Stochastic Queue Median (SQM) model. SQM assumes that if calls are waiting in 
queue to be served, the dispatcher, when picking the next call to be serviced, will 
consider all calls to be of the same importance and will choose the next call by a rule 
that depends only upon the relative position of the calls in queue. For example, SQM 
allows for queueing disciplines like First-Come-First-Served (FCFS), Last-Come-First- 
Served (LCFS) and Service-In-Random-Order (SIRO), etc. 

This type of priority rule is clearly inappropriate in many contexts. For example, 
in urban emergency services, calls that involve danger to human life (e.g., a major 
accident in the case of an ambulance system or a “violent crime in progress” in the 
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case of police patrols) deserve and receive higher priority over calls for more routine 
incidents. A realistic model for such contexts should allow for prioritization of calls 
for service. 

Our aim in this paper is to formulate and provide solution techniques for a single 
server queueing-location model which allow calls to be selected from an arbitrary 
number of priority classes. The exact number of priority classes which should be 
employed depends on the particular application at hand and requires study in its own 
right. 

The reason for restricting our attention to a single server model is that, in the course 
of our analysis, we shall need a closed-form expression from queueing theory for the 
average waiting time in queue for calls of various priorities. While such an expression 
does exist for an M/G/1 queueing system, as of today, none exists for an M/G/m 
system when m (the number of servers) is greater than one. Two attempts to solve 
multiple server queueing-location models through an approximation technique are the 
papers by Berman and Larson [2] and Larson [9]. 

We refer to the model developed in this paper as a K-Priority Queueing-Location 
(K-PQL) model. The principal results presented in this paper can be summarized as 
follows: 


1. The optimal K-PQL model location is usually different from that obtained by 
grouping calls from all priorities into a single category and using the SQM model 
developed in [3]. 

2. For a tree network we establish properties of the objective function that enable 
one to trim off portions of the tree where the optimal location cannot lie. 

3. We characterize the set of optimal locations for extreme values of the arrival 
rate. For a tree network these exist at unique points. For a cyclic network the set has 
finite cardinality. Elements of the set can be efficiently identified in both cases. 


FORMULATION 
Given an undirected network G = (N,A), with N being the node set, | N | =n, A 
being the arc set, | A | = a, assume that travel between two points x and y on G 


occurs on a shortest length path connecting these two points and denote the length of 
the path by d(z,y). 

Calls arrive solely at the nodes of the network. They are divided into K priority 
classes, indexed by k, k = 1,2, ... ,K; the lower the value of k, the higher the 
priority of the call. Arrivals of calls of priority & constitute arrivals from a time 
homogeneous Poisson stream, with average arrival rate A,. Given that an incoming 
call is of priority k, the a priori probability that the call originated from node i is fix. 
We also define A = D*_,A, and F, = Aj,/A. 

The service system operates as follows. The server is located at a point on the 
network, called the “home location,” when not busy servicing a call. After servic- 
ing a call the server always returns to the home location, before answering the next 
call. If a call enters the system, and the server is busy serving another call, the 
arriving call is entered into a queue. The queue is depleted by employing a Head- 
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of-the-Line (HOL) non-preemptive priority queueing discipling, FCFS within each 
priority class. 

The system time of a call can be divided into two portions: (a) queueing delay; and 
(b) service time (travel time to call, on-scene service time at scene of call, travel time 
from call, and off-scene service time at home location). 

The server, when travelling to a priority & call, travels with constant speed v,. 
Furthermore, when he travels from the scene of the priority k call to his home location, 
he travels with speed v,/(® — 1), B > 1, where B is a prespecified constant. 

The first and second moments of nontravel time related service time (on-scene plus 
off-scene) for a priority k call from node i, are denoted by finite constants W,, and 
W3,, respectively. In particular, we note that these are assumed independent of the 
home location of the server. 

To optimize for the “best” home location of the server we have to decide upon an 
objective function. In an urban emergency service system and for a priority & call, 
quality of service if often measured with reference to Response Time, TR,. When the 
home location is at x on the network, 


TR(x) = Qi) + 400, (1) 


where 
Q,{x) = queueing delay of a priority & call, when 
the server’s home location is at x, and 
t,{x) = travel time to the priority k call, when the 
server’s home location is at x. 


TR,(x), O,(x), and 4,(x) are random variables, as a result of the spatial and temporal 
uncertainty associated with the calls for service. Denoting the average of a random 
variable by y by y, and attaching an importance factor 1,(7™, = 0) to priority & calls, 
the weighted average response time to a call, when the server’s home location is at 
x, is given by: 


K 
TR(x) = >) mu TR,(2). (2) 


k=1 


We wish to find a point x which minimizes TR(x). 


ANALYSIS 


To simplify presentation, the following notation shall be used. A variable A,(x) 
indicates the variable A depends on (i) the priority class k, and (ii) the home location 
x of the server. A variable B(x) indicates that variable B depends on the home location 
x of the server. 

Denoting the expected travel time to a call, expected service time and second moment 
of service time by the symbols 1, §, and S?, respectively, we have in terms of variables 
previously introduced: 
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ix) = Dd fa dx.d/v, (3) 
iEN 
Sx) = D> fa (Wa + Bdlx,d/vi, (4) 
iEN 
Sx) = D falWk + 2WaBdx,d/v, 
iEN 
+ (Bd(x,i)/v,)"1, (5) 
— K — 
Sx) = DS FSX). (6) 


k=1 


The assumption that the server always returns to his home location before answering 
successive calls leads to independent and identically distributed (i.i.d.) service times. 
The system operates as an M/G/I queueing system with K priority classes, non- 
preemptive priorities, FCFS within each priority class. 

The average queueing delay for priority k calls, denoted by symbol O, is (for 
example, see Kleinrock [8]): 


Qua) = 4 >: nie 500 | 


if 1 — S*_, ASfx) > 0, (7) 


L + %, otherwise. 


Using (1)+{7) we can therefore write an expression for the weighted average response 
time, TR(x). a 

In order to search efficiently for the minimum of TR(x), we identify regions of the 
network where TR(x) takes on a unique functional form. We briefly outline here the 
approach used, as it is virtually identical to the one in [3] and Chiu [4]. The weighted 
travel distance function is piecewise linear and concave over any arc of a network. 
We will define the region of an arc over which this weighted distance function is linear 
as a primary region. Further, various associated functions (when finite) are differen- 
tiable over the interior of a primary region. _ 

The steps in identifying a minimum of TR(x) on the network are: 


(1) Identify the primary regions on the network. 

(2) Eliminate primary regions (for a tree network) where the minimum of TR(x) 
cannot lie. 

(3) Find a local minimum for each of the remaining primary regions. A descent 
algorithm (for example see Shapiro [10]) can be used for this purpose, after the feasible 
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FIG. 1. IHlustration for proof of Theorem 3: (a) Graph G, (b) Subgraph G'. 
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portion of the primary region has been identified. The feasible portion of a primary 


region is the section on which TR(x) is finite. 
(4) Compare these local minima to obtain the home location x on G. 


SIMPLIFICATION FOR TREES 


In contrast to the case of cyclic networks, we can establish properties of the objective 
function on trees which allow us to eliminate primary regions or portions thereof, 


where the minimum of TR(x) cannot lie. This is the purpose of this section. 


Theorem 1. Q,(x) is strongly quasiconvex (when finite) on any path of a tree network, 


forA > 0, Vk = 1,2,... ,K. 


Proof. See Batta [1]. 


Theorem 2. Sf ,712;(x) is convex on any path of a tree. 


Proof. Convexity of t,(x) on any path of a tree is well known, see Goldman [6]. 


The result follows from the fact that 7, = 0, Wk = 1,2,... ,K. 
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We use the results stated in Theorems 1! and 2 to eliminate portions of a tree network 
where the global minimum for TR(x) cannot lie. 


Theorem 3. Let x,,x2, . . . 4x be the unique points on the tree network G which 
Tespectively minimize 2100), Q(x), .. . , Ox(x). Let g be any point which minimizes 
ZK_ Tt, (x). Finally, let G’ be the smallest connected subgraph of G which contains 
all the points, x,, x2, ... Xx and g. Then the minimizer of TR(x) lies on subgraph 
G'. 


Proof. The proof is by contradiction. Suppose this were not the case. Then for a 
tree network G (Fig. 1a) and associated subgraph G’ (Fig. 1b) a minimizer of TR(x) 
lies at a point x* € G’, as shown in Figure 1a. 

By moving the location from x* to the nearest point on subgraph G’ (node q in Fig. 
la), O:(x), Qo(x), . . . , Og(x) all strictly decrease, by Theorem 1. Also DX. tty (x) 
cannot increase, by Theorem 2. From (1), (2) and from the fact that 7, = OWk = 
1,2, ... ,K, it follows that TR(x) must strictly decrease through this movement. The 
Theorem follows. | 


To complete our discussion on elimination of primary regions of a tree network G 
where the minimum of 7R(x) cannot lie, we state algorithms for finding: 

(a) A minimizer of Df. tral, (x). _ 

(b) Minimizers of Q,(x), Q2(x), . - - » Ox(x). 

The minimizers found in (a) and (b) can then be used to form the subgraph G’ of 
G to which we can restrict our search. 


Algorithm for Finding Minimizer of Di. ,714t.(x) 


A nodal minimizer of =f_ ,74f, (x) can be found by using Goldman’s algorithm for 
finding the one median on a tree, see Goldman [6], with the weight of node i being 


_vx« 
wie DeeiT i Me 


Algorithm for Finding Minimizer of Q,(x), kK =1,2,...,K 


By Theorem 1, O,(x) is strongly quasiconvex (when finite) on any path of a tree, 
Vk = 1,2,... ,X. An algorithm very similar to that used in [5] (for finding the 
SQM on a tree) can be developed. We omit the details. 


SENSITIVITY ANALYSIS WITH RESPECT TO ARRIVAL RATE 


In this section, we address the question: How does the optimal location of the facility 
change as the arrival rate of calls is increased from 0 to the maximum possible arrival 
tate, Amax (Amax is the minimum value of \ which renders every point of the network 
infeasible as a home location; specifically, if node j is a minimizer of =E_ ALS, (x), 
then Amax = W/(SK_ iF, S,(j)). We formalize our observations in Theorem 4. 
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Theorem 4. Define sets S, T, S, and 7, as follows: 


S= {. | x€EG, Se Td lX) — Dr. mde), VyEe of. 


x|x€S, M(x) < eel 


ts 
Hi 


x|x€G, DK NSX) S = MS), Wy of, 


| 
ee cc: 


T; x | x€ Se - =-_ =... = = 
E = Mineo "Fe 54) 


Sy) 


a, Vv ES. 
|: = Amax nr Fe 5409 


IA 


Then, the following hold: 

(a) for every Xo & T and yo € T, one can find e, > 0 such that TR(xo) > TR(yo) 
whenever O = A S €), ie — 

(b) for every Xo€ T, and yo € T,, one can find e, > O such that TR(xo) > TR(yo) 
whenever Amax — €2 SA < Amax, and 

(c) the sets T and T, have finite cardinality. When G is a tree T and 7, are singleton 
sets. 


Proof. See Appendix. e 


To present intuitively the contents of Theorem 4 we distinguish two kinds of Hakimi 
medians. A point is a type ] Hakimi median, denoted by HM], if it minimizes 
DK_ yet, (x) on G. Similarly, a point is a type II Hakimi median, denoted by HM2, 
if it minimizes DK_ ,A,S, (x) on G. HM1 is the set of points which minimizes the 
weighted (by priority class) average travel time to a call, and HM2 is the set of points 
which minimizes the system average service time (or system average travel time) to 
a call. 

If HMI is unique, then from [7] is exists at a node j of the network (the fact that 
at,’S are not nodal weights does not change this observation). In this case, Theorem 
4 asserts that node j is the optimal location for A close to 0. If HM1 is not unique, 
then the optimal location is a point in the set of points HMI, which minimizes Stx). 

If HM2 is unique, then from [7] is exists at a node k of the network. In this case, 
Theorem 4 asserts that node & is the optimal location for \ close to Ajax. If HM2 is 
not unique, then the optimal location is a point in the set of points HM2, which 
minimizes a function of S*(x)/D, where D depends only on the calls of the first K — 
1 priorities. 

Theorem 4 also asserts that the optimal location is unique when d is close to either 
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O or Amax for the case of a tree network, and the set of optimal locations is finite for 
a cyclic network. 

To find the elements of 7 and 7, on a tree network we can exploit the strict convexity 
of S7(x) and 


$*(x) 
E — max Deo Fs 54 | 


on any path of a tree, to develop algorithms similar to that of finding the minimizer 
of Q,(x) on a tree. 

To find elements of T and 7, respectively on a cyclic network we can exploit the 
strict convexity of the above mentioned functions on a primary region. 

For intermediate values of \, the average queueing delay of a call and the average 
travel time to a call are of comparable magnitude. The second moment of service time 
plays a central role in determination of the optimal location. Our computational ex- 
perience with the model indicates similar behavior to that observed in [3] and [5]. 


NUMERICAL EXAMPLE 


In this section we illustrate the results of this paper through a simple numerical 
example. 

We use the problem instance described in Figure 2. The first column of Table I 
displays the optimal location when ) varies from 0.045 to 0.855 under the 2-priority 
queueing-location (2-PQL) model. The second column of Table 1 displays the corre- 
sponding optimal location obtained by merging the calls from priorities 1 and 2 into 
one category, and using the SQM model in [3]. The third column of Table I gives the 
corresponding optimal location when ignoring queueing and using the model in [7], 
which we refer to as the Hakimi model. 

Table I displays the objective function value obtained under the 2-PQL model when 
the facility is located at (i) the optimal 2-PQL model location, (ii) the optimal SQM 
model location, and (iii) the optimal Hakimi model location. Define: 


A(x) = tili(x) + Tratd(x), 


B(x) = F\S\(x) + F2S,(x). 


A. We note that A(x), the weighted sum of average travel times to calls of priorities 
1 and 2, is uniquely minimized at node 1. 

B. B(x), the system average service time to a call (of priority | or 2), is mini- 
mized_at any point on arc (3,5). 

Cc. S*(x), the second moment of the service time of all calls (of priority 1 or 2) is 
minimized uniquely at a point 2.5 units from node 3, on arc (3,5). Here we have an 
unusual situation: the point which minimizes the second moment of service time also 
minimizes average service time. 


D. Amex = WFS\(x) + F2S2(x)) 


x is on are (3,5) 
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FIG. 2. Numerical example. 

Number(s) next to arc(s) represent arc lengths. B = 2; travel to a call is with same speed as 
travel from a call. 7, = 2; m2 = 1; importance factors for priority 1 and 2 calls respectively. 
P, = [1/2,1/8,1/8,1/8,1/8]; i th entry is probability that the next priority | call comes from 
node i. P2 = [1/8,1/8,1/8,1/8,1/2]; i th entry is probability that the next priority 2 call comes 
from node i. v; = v2 = 10; the travel speeds for priority 1 and 2 calls respectively. W, = 
(0.1,0.1,0.1,0.4,0.1]; i th entry is mean non-travel time related service time for priority 1 calls 
from node i. W, = [0.1,0.1,0.1,0.1,Q.1); i th entry is mean non-travel time related service 
time for priority 2 calls from node i. W? = [0.01,0.01,0.01,0.01,0.01]; i th entry is second 
moment of non-travel time related service time for priority 1 calls from node i. W3 = 
(0.01,0.01,0.01 ,0.01,0.01); i th entry is second moment of non-travel time related service time 
for priority 2 calls from node i. F, = 1/3, F, = 2/3; fraction of calls of priority 1 and 2 
respectively. 


Further define: = = 
C(x) = S? (xy (1 — Xanax FrSi(2)). 


Then C(x), when x is on arc (3,5), is minimized uniquely at a point 1.77 units from 
node 3, on arc (3,5). From the assertion in Theorem 4(b), the optimal 2-PQL location 
is at the point which minimizes C(x) on arc (3,5), when A is sufficiently close to Amax- 
Now we are in a position to discuss the results in Tables I and II. 
In our example the minimizer of S?(x) also minimizes B(x), therefore we are assured 


TABLE I. Optimal location under (a) Hakimi Model, (b) SQM Model, 
and (c) 2-PQL Model. 


SQM Model (priorities 


r 2-PQL Model grouped together) Hakimi Model 

0.045 node 1 arc (3,5), 2.5 from any point on arc 
node 3 (3,5) 

0.180 node 1] arc (3,5), 2.5 from any point on arc 
node 3 (3,5) 

0.360 arc (3,5), 0.995 arc (3,5), 2.5 from any point on arc 
from node 3 node 3 (3,5) 

0.540 arc (3,5), 1.451 arc (3,5), 2.5 from any point on arc 
from node 3 node 3 (3,5) 

0.720 are (3,5), 1.598 arc (3,5), 2.5 from any point on arc 


from node 3 node 3 (3,5) 
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TABLE II. Objective function value when using 2-PQL Model assumptions and locating 
facility as in Table I. 


Hakimi Model 
SQM Model (priorities (assuming location 

nN 2-PQL Model grouped together) is at node 3) 
0.045 1.461 1.760 1.589 
0.180 1.973 2.148 2.031 
0.360 2.813 2.882 2.842 
0.540 4.117 4.186 4.242 
0.720 7.827 7.955 8.199 


that for all feasible values of \ (0 < \ < Amax), the optimal location under SQM will 
be a distance 2.5 from node 3, on arc (3,5). 

For the deterministic case the optimal location is any point on arc (3,5); we have 
chosen node 3 arbitrarily. 

For low values of A, the 2-PQL model favors node | as an optimal location primarily 
because of the concern reflected for priority | calls (7, = 2, m7) = 1). 

For intermediate values of A, queueing delays and travel times to calls are of 
comparable magnitude: hence the optima! location for the 2-PQL model shifts towards 
the minimizer of the second moment of service time. We note that the higher importance 
attached to priority 1 calls tends to keep the location toward node 1, hence the optimal 
location is closer to node | than under the SQM model. 

For large values of }, we appeal to Theorem 2(b) to yield the optimal location of 
the facility under the 2-PQL model. 


CONCLUDING REMARK 


We end this paper by commenting on properties of the algorithms developed therein. 

TR(x) is finite if. Q,(x), Q2(x), . . . , Qx(x) are finite. From Lemma | (in Appendix) 
we know that the Q,(x)’s are strongly quasiconvex (when finite) on a primary region 
of the network. Even though 7R(x) is a positively weighted sum of the O,(x)’s this 
does not necessarily imply any convexity property of TR(x) (when finite). 

The strong quasiconvexity (when finite) of TR(x) persists in numerical examples 
because (see proof of Lemma 1 in Appendix) the only case when Q;(x) may not be 
strictly convex (when finite) in a primary region is when its denominator is strictly 
convex; this happens when the mean service times with respect to calls of the first k 
— 1 and k priority classes are uniquely minimized at the same end of the primary 
region. 


The authors would like to thank the anonymous referees for their helpful comments on an 
earlier version of this paper. This research was supported, in part, by the National Science 
Foundation, viz-a-vis Grant No. ECS-8204318; the support is gratefully acknowledged. 


APPENDIX 


This Appendix contains proofs of some of the theorems and lemmas in this paper. 
In order to prove Theorem | we need to prove a Lemma on strong quasiconvexity of 
Q,(x) on a primary region of a network. 
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x primary region 


2 (>) 


FIG. 3. Illustration used in proof of Lemma 1. 


Lemma 1. Q,(x) is strongly quasiconvex (when finite) on a primary region of a (cyclic 
or tree) network, for’ > 0, Vk = 1,2,... ,K. 


Proof. Consider a primary region 1-2 on arc (a,b) of the network, as shown in 
Figure 3. > 

We wish to prove strong quasiconvexity of Q,(x,), A > 0, fork = 1,2,...,K, 
when the function is finite. 
_ For notational convenience we shall suppress the parameter x in the formulae of 
Sx), S2(x), (x) and $?(x), labeling them as S,, SZ, ,, and S?, respectively, whenever 
the context leaves no ambiguity. In addition, we shall denote first and second derivatives 
of a function f(x) by f'(x) and f"(x); further suppressing x, we write these as f’ and 
f’, Tespectively. 

We divide the analysis into two cases: 


Case 1; k = 1. In this case, the denominator of O1(x) is a linear function of x. It 
follows from [3] that Q,(x) is strictly convex in a primary region. 


Case ll: k > 1. 
= (n/2)S2 
Q(x) = k-1 


[25] [+- 25] : 


Let N(x) = (N/2) S = a + bx + cx’, and 


D(x) 


k-1 k 
E cig 5; | E = > 5 | = (a, + b,x)(a. + box) 
j=l 


j=! 


=dt+ext fr, 


where a,b,c,d,e, and f are constants that can be evaluated for a given primary region. 
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Since \ > 0, it follows from [3] that: N(x) >OWx ERD b* - 4ac<0> 
(N’)? — 2N'N <0. We note that N, N’, and N” are defined for x = 0. 

Also, D(x) has real roots, because D(x) = (a, + b,x) (@2 + box). 

Direct differentiation of Q,(x) yields: Qy(x) = { D}°{N"D? -2 D'DN' 
+ 2(D'.N — ND"D}. 

Since Q,(x) is assumed finite, it takes on the sign of: 


T(x) = N’D? — 2 D’DN’ + 2(D’)? N — ND’D. 


We distinguish 3 cases. 


Subcase (a): D” = 0. in this case D(x) is linear in x and the arguments of Case I 
apply. 


Subcase (b): D’ < 0. In this case, we can write: 
T(x) = N"D? — 2D'DN' + 2(D'YN — ND"D 
2 N"D*? — 2D'DN' + 2(D'YPN, 
where the inequality follows from N,D > 0 and D” < 0. The proof that 7(x) = 0 V 


x € R follows directly from Case I, because D” has been effectively removed from 
the function, and the proof is now identical to that presented in [3]. 


Subcase (c): D" > 0. We distinguish two cases, when the roots of D(x) are distinct, 
and when they are not distinct. Throughout this discussion we shall temporarily assume 
that Q,(x) is defined by (8) even if D(x) < 0. 


Sub-subcase (i): roots of D(x) are distinct (Fig. 4). Let the roots of D(x) be y, and 
Y2, respectively. In the interval (y,,y2), D(x) is negative. Therefore, Q(x) takes on 
the sign of —7/x). 


—T(x) = —N"D? + 2D'DN' — 2(D'YN + ND'"D 
< —N"D? + 2D'DN' — 2(D'yN, 
where, the last inequality follows from N,D" > 0 and D < 0. 
Let G(x) = —N'D? + 2D'DN' — 2(D'YN 
= —N"D? — [2D'DN' + 2(D')N]. 
It follows that G(x) < 0 V x € R, from the arguments in subcase (b). Therefore Ox) 
is strictly concave in x, when x € (y1,y2). _ 

Qi{x) = 0, yields a quadratic in x (the cubic terms cancel out). Since Q,(x) > —®© 
when x— y, from above, and Q,(x) > —® when x— y; from below, strict concavity 
of Q,(x) in the interval (y,,y2) implies that exactly one of the solutions to Qj(x) = 0 
lies in the interval (y,,y2). Therefore, exactly one solution to Qj(x) = O lies in R — 


(y1,¥2), where R = (—®, +), This establishes strong quasiconvexity of O.(x) when 
D(x) > 0. 
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Illustration for proof of Lemma 1, Subcase (c): Sub-subcase (i). 


FIG. 4. 
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Illustration for proof of Lemma 1, Subcase (c): Sub-subcase (ii). 


FIG. 5. 
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Sub-subcase (ii): roots of D(x) are not distinct (Fig. 5). Let the double root of D(x) 
be x*. Then exactly one solution to Q;(x) = 0 lies in the interval (—©,x*], and one 
solution in the interval [x*,%). 

Strong quasiconvexity of Q,(x) in the intervals (— ©,x*) and (x*,®) follows from 
the argument in Sub-subcase (i). 

The lemma follows. a 


Proof of Theorem 4. (a) The proof is by contradiction. Suppose 3 xo & T and yo 
€ T, such that 7R(xo) < TR(yo), for 0 < X_S &,, for every €, > 0. 

Initially pick €; = Amax/2 > O. Then 7TR(x) is finite for O =< A < e€,, and its 
expression becomes: 


TR(x) = (11) 


Nin 


+ 1,(x). 


$*(x) 
= [ - ry F540) | [! - E50) 


j=l 


We consider two cases. 
Xo € S — T. Here xo,yo € S, and Dko yyy (Xo) = Thome tyo)- 
Upon substituting a = S?(x)/2, by = Diz} F)S{xo), 


= Dh FS A%0), 
d = S%(yo)/2, ey = Zkz}F;S; (yo) and 
& = Zi_F)S(yo), we get: 


(note that A, by, Cyd, 2 >OVk= 1,2, areca »K) TR(Xo) = TR(yo) 


= ADf-1% 
& — d) + N-ale, + fi) + dy + cd) + Mla he — al (12) 
(1 — Ad, J{1 — Ac,J{1 — Ae]1 — Agi] 
= AMA) < 0. 


We note thata > d, 1 — AB, > 0, 1 — Ae > O, 1 -— Ae, > O, 1 — AQ, > O and 
uw, > 0, fork = 1,2, ... ,K. Therefore (12) leads to a contradiction when e, is 
chosen small enough, 0 < €, S Amax/2 andO <A < e&). 

Xo ES — T. In this case, xo € S, and Dk, mitk(Xo) > Thy Mak; (Vo). 

Using notation introduced earlier, we can write: 


K K 
TR(xo) — TR(yo) = DS tddxo) — DS tudlyo) + MOA) < 0. (13) 


k=1 k=] 


Let p* = ming<,<e, (A), and s* = min(—1,p*). 
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K K 
~ atruti(X0) — ard 
ieee na [> iti{Xo) 2 uo | Res > 0, 


2s* 7 2 


Then (12) is a contradiction for 0 < A < &. 
(b) The proof is once again by contradiction. Suppose 3x» & T, and yo € T,, such 
that TR(xo) < TR(yo), for Xinax — €2 = A_< Amar for every €2 > 0. 


Initially pick €2 = Amax/2 > 0. Then TR(x) i is as in (11). We consider two cases. 
x ES, — Ty. 
Let 
K-1 
g(A) = 2 Tr 


E -— d) + N-ale, + f,) + dh + c)) + Mlaeh — Pe 
{] — Ab,)[1 — Ac,J{1 — Aexjl — Ag] 


L 


K K 
max dg(d) + DS) mdxo) — DS tete(yo), and 


O=A<Amax k=l k=! 
L* = max(1,L). 
Then TR(Xo) - TR(yo) < 0 implies that 


Aa Ad 
ae A li Shed eS <a ao 


Since x» € S, — T,, we have cx = gx, therefore (14) reduces to: 


ATK a d 
L* + —— | -————_—— - ———] < 0. 
1 — Ack Ee — dbx] {1 - <aal ~ 


Since x) © S; — 7, a@/(1 — AmaxPx) > d/l — Amaxx), Which implies a/(1 — Abx) 
> dil — dex), OS A < Amax. Let T* = a/(1 — AmaxDx /2) — d/i(L — Amaxex/2) > 
0. Then (15) reduces to: 


d 
pig a 2g. (16) 
1 - ACK 


Set ©) = Amax — mMax[Amax/2, L*/(L*cx — tx T*)] < Amax, Where the inequality 
follows because 1x,7* > 0. We note that 2 > 0. For Amax — €2 S A < Amax, (16) 
yields a contradiction. 
Xo ZS, — Ty. Pick 82 = Weg < 1 lex = Amar Then TR(x0) becomes + © for 8, 
=D < Amax, Whereas TR(yo) is finite. Therefore TR(xo) _ TR(yo) < 0 would yield a 
contradiction for Xmax — &2 <A < Amax, Where €2 = Amax — 42. 
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(c) The proof of this follows from the following facts: 

(i) S*(x) is strictly convex on any path of a tree, and on a primary region of a cyclic 
network. 

(ii) SGV. — Amax ZZ} F,S;(x)) is strictly convex (when finite) on any path of a 
tree, and on a primary region of a cyclic network. 

(iii) For a tree network if x,y € S, then every point on the path p/x,y/ joining x 
and y belongs to S. The same applies for S,. 

(iv) § and S, consist of a finite collection of connected intervals on a cyclic network. 

The theorem follows. a 
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